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ABSTRACT 



A flow visualization study of natural convection in a 
liquid filled rectangular enclosure with a small heater pro 
truding into it from one vertical wall was conducted. The 
top and bottom horizontal surfaces of the enclosure were 
heat exchangers whose temperatures could be varied independ 
ently. The fluid in the enclosure was water. The Baker 
electrochemical technique which utilizes a pH indicator was 
used for flow visualization. Photographs were taken of the 
two dimensional flow for three different locations of the 
heater on the vertical wall. Nusselt numbers for each 
heater location and for a range of Rayleigh numbers were 
also determined. Visual results indicate the presence of 
a dual-cell configuration within the enclosure. There is 
a buoyancy driven upper cell and a shear driven lower cell. 
The flow of the upper cell follows the geometry of the 
enclosure. Though the uncertainty inherent in the experi- 
ment was high data calculations suggest a trend that as the 
heater is raised within the enclosure the Nusselt number 
decreases . 
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NOMENCLATURE 



Symbol 


Description 


Units 


A 


Area 


.2 .^2 
in or ft 


g 


Acceleration o£ gravity 


2 

Ft/sec 


Gr 


Grashof number 


Dimensionless 


h 


Convection heat transfer 
coefficient 


Btu/hr-ft“-°F 


k 


Thermal conductivity 


Btu/hr- ft-°F 




Fluid thermal conductivity 


Btu/hr-ft-°F 


^g 


Thermal conductivity of plexiglas 


Btu/hr-ft-°F 




Thermal conductivity of foam 
rubber insulation 


Btu/hr-ft-°F 


>'p 


Thermal conductivity of phenolic 


Btu/hr-ft-°F 


>=ss 


Thermal conductivity of 304 
stainless steel 


Btu/hr-ft-°F 


L 


Characteristic length 


Inches 


Nu 


Nusselt number 


Dimensionless 


Pr 


Prandtl number 


Dimensionless 


^cond 


Heat loss via conduction through 
vertical walls 


Btu/hr 


Q 

^conv 


Energy convected into fluid 
from heater 


Btu/hr 


Qin 


Energy into cartridge heater 


Btu/hr 


R 


Resistance of precision resistor 


Ohms 


Ra 


Rayleigh number 


Dimens ionless 
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Equivalent thermal resistance to 
conduction through viewing wall 


Hr-°F/Btu 




Equivalent thermal resistance to 
conduction through wall opposite 
viewing wall 


Hr-°F/Btu 




Equivalent thermal resistance to 
conduction through wall holding 
heater 


Hr-°F/Btu 


^cond 


Total thermal resistance to 
conduction with respect to heat 
loss through vertical walls 


Hr-°F/Btu 




Individual thermal resistance 
to conduction Cj=l>7) 


Hr-°F/Btu 




Ambient temperature 


°F 




Average temperature of heat 
exchangers 


°F 


Th 


Cartridge heater temperature 


°F 


Ti 


Thermocouple temperatures Ci=l,8) 


°F 


Ts 


Average temperature of heater 


°F 




Voltage across heater 


Volts 


^2 


Voltage across precision resistor 


Volts 


W 


Uncertainty of variables 


Various 


a 


Thermal diffusivity 


2 

Ft /sec 


B 


Coefficient of expansion 


1/°F 


V 


Kinematic viscosity 


Ft^/sec 
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I . BACKGROUND 



A. INTRODUCTION: STATEMENT OF PROBLEM 

The rapid advancements in the field of miniaturization 
of electronic components is well-known. Technology has 
brought us from the large-scale integration (LSI) device to 
the very large-scale integration (VLSI) device in just a feiv 
short years. Since the first 1 kilobit RAM semiconductor 
device was introduced in 1970 the number of memory cells on 
a single device has grown to 64K while the size of the device 
has decreased. One such device, the 64K RAM Micron chip pro- 
duced by Tiny Micron Technology, Inc. is only 22,000 square 
mils, or approximately the size of one of these printed 
letters [1] . Yet microminiaturization advancements are far 
from over -- 256K and greater RAM chips are on the immediate 
horizon . 

As in most rapidly developing technologies there are var- 
ious problems associated with the advancements that have been 
made in the miniaturization of electronic devices. One of 
the primary problems that has evolved is that of suitably 
cooling components so that they will perform reliably over 
long periods of time. For though the size of the individual 
component has shrunk considerably its junction temperature 
has remained fairly constant. It follows that the decrease 
in device size has been kin to an increase in heat flux. 
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To get an idea o£ the magnitude o£ this decreased size/ 

increased heat £lux relationship one only has to heed the 

words o£ Chu [2]. He asserts that the heat £lux o£ a micro- 

circuit chip is but two orders o£ magnitude less than the 

5 7 

heat £lux o£ the sun (10 vice 10 Watts per meter squared) . 
Further, the microcircuit chip must dissipate its energy at 
strict maximum junction temperatures o£ between 125 to 150 
degrees C. The integration o£ a multitude o£ devices on a 
single printed circuit board leads to an increase in the 
temperatures in which the device operates. The high heat 
£lux o£ the microcircuit chip coupled with the increase in 
environmental temperature very o£ten leads to the problem 
o£ component £ailure resulting £rom inadequate cooling. It 
is this problem that has recently prompted extensive research 
into the techniques o£ cooling microelectronic devices. 

B. DESIGN CONSIDERATIONS; PROPOSED SOLUTIONS 

Many papers have been published outlining alternative 
methods o£ thermally controlling microelectronic devices. 
Hannemann [3] describes the thermal considerations that must 
be taken into account at each level (subassembly, assembly, 
and system levels) o£ the device design process. Kraus, et 
al [4] assert the need £or thermal analysis even a£ter 
design. They propose a strategy that yields conservative 
estimates o£ the heat £lux associated with electronic cir- 
cuitry. Further, Bar-Cohen, et al [5] indicate that numerical 
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analysis of the heat transfer problem should be relied upon 
to cope with the analytically intractable nature of even 
standard packaging configurations. 

Once the initial design step has been completed a reli- 
able cooling scheme must be employed. A number of such 
schemes including conduction, forced air cooling, and direct 
liquid cooling have been proposed by many authors. Bergles , 
et al [6] contend that of the many techniques available, 
direct liquid cooling may be the most desirable mode of heat 
transfer when considerations such as weight, size, or cost 
are taken equally into account with the heat transfer coef- 
ficient. Simons and Moran [7] support this view by stating 
that liquid coolants offer the capability to achieve heat 
transfer coefficients one to two orders of magnitude higher 
than that of air. 

Many authors have offered alternative methods of utiliz- 
ing immersion cooling. One such scheme is proposed by Kraus 
and Bar-Cohen [8] who analyze submerged condensers. They 
describe a model whereby the heat from the immersed elec- 
tronic devices is removed via bubble pumped augmented con- 
vection interacting at the submerged condenser surface. 

This study focuses on direct liquid cooling by natural 
convection. More specifically, since most microelectronic 
devices are seated on printed circuit boards packaged very 
near each other, investigation is limited to natural convec- 
tion in a small liquid filled rectangular enclosure. One 
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advantage of natural convection cooling is that there are no 
pumps, pipes, or any other equipment within the enclosure; 
the manufacture of a relatively small module may be assured. 

C. IMMERSION COOLING: ANALYTICAL STUDIES 

Numerous analytical studies of natural convection in rec- 
tangular enclosures have been made since 1966 when Wilkes 
and Churchill [9] initiated finite difference analysis for 
such configurations. This work was followed by many others 
including Szekely and Todd [10] who, in 1970, presented a 
finite difference analysis of natural convection in a rec- 
tangular cavity in which two opposing vertical walls provide 
the temperature gradient. Comparing their numerical results 
with experiment they concluded that within the Rayleigh num- 
ber range of sixteen hundred to s ixty- thousand the thermal 
response of a system is approximately independent of the 
Rayleigh number. They further noted a trend for secondary 
cells to develop once the Rayleigh number exceeded three 
hundred thousand. 

In 1975 Chu, et al [11] studied, via finite difference 
methods and experiment, the natural convection in a rectan- 
gular enclosure in which the heat source was a long, hori- 
zontal strip located in one vertical wall. Focusing on 
heater location, heater size, and enclosure size a conclu- 
sion was reached that as the Rayleigh number increases the 
heater location for maximum circulation shifts downward but 
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that heater size seems to have little effect on circulation. 
They contend that secondary cells develop at the upper sur- 
face as the height of the enclosure is increased. 

Markatos and Pericleous [12] have recently published a 
study of buoyancy-driven flow in a square cavity with dif- 
ferentially heated sides. Secondary cells were also observed 
in their study. 

D. IMMERSION COOLING: EXPERIMENTAL STUDIES 

Megerlin and Vingerhoet [13] collected natural convection 
data using fluorocarbon dielectric fluids in an enclosure 
modeled after a microelectronic module. One significant con- 
clusion reached was that cold plate orientation has a signi- 
ficant effect on the heat transfer --a fifty percent increase 
is attained using a horizontal vice vertical cold plate. An 
experimental investigation using an undivided and partially 
divided water filled two-dimensional rectangular enclosure 
was performed by Nansteel and Greif [14] . 

In 1972 Baker [15] published the results of his experi- 
mental work in which he used resistors immersed in Freon. 

He concluded that the heat dissipation obtained with Freon 
is at least four times greater than that obtained using the 
same assembly with stagnant air. Lin and Akins [16] investi- 
gated natural convection heat transfer in a fluid filled cube 
where all six wall temperatures were subjected to a step 

change in temperature. Their conclusion was that a variety of 
circulation patterns evolve depending on the size of the cube. 
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In 1980 Ostrach, et at [17] published the results o£ an 
experiment that used small particles introduced in silicone 
oils that filled various sized rectangular enclosures with 
height to length ratios less than unity. Secondary cells 
were evident in their study. Finally, an experiment was per- 
formed by Bohn and Kirkpatrick [18] which studied high Rayleigh 
number natural convection in a cublical enclosure filled with 
deionized water using various combinations of heated walls. 

They concluded that heating from below along with simultaneous 
vertical wall heating and cooling enhances heat transfer. 

The above discussion is not intended to be a definitive 
synopsis of all works in the design considerations of micro- 
electronic devices, nor in the numerical analysis and experi- 
ments conducted on natural convection in enclosures. References 
cited are merely intended to be representative of the works 
that have been done in these areas. They serve to provide an 
excellent point of departure for the remainder of this study. 

For a more extensive review of the works done concerning natu- 
ral convection in rectangular enclosures the reader is referred 
to Catton [19 ] . 

E. OBJECTIVES 

The objectives of this study are; 

1. to design and build an apparatus that approximates a 
liquid filled microelectronic package with appropriate printed 
circuit board spacing that is cooled by heat exchangers acting 
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as the top and bottom horizontal surfaces. The heat source 
is to be a rectangular protrusion modeled geometrically after 
a 40-pin dual-in-line package (DIP) that is attached to and 
horizontally traverses the length of one vertical wall. 

2. to observe the steady state natural convection flow 
in the apparatus using a variety of heater and heat sink 
temperatures as well as various locations of the heating 
element . 

3. to analyze data obtained and determine Nusselt and 
Rayleigh numbers for each apparatus configuration. 

This study is intended to serve as a basis for future 
heat transfer experiments in natural convection using vari- 
ous enclosure sizes and multiple heated protrusions. 
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II. EXPERIMENT 



A. GENERAL DESIGN CONSIDERATIONS 

The flow visualization apparatus (Fig- 1) consists of a 
small liquid filled rectangular enclosure with a heater 
modeled geometrically after a 40-pin DIP protruding into 
the enclosure from one of the vertical walls. This con- 
figuration is meant to simulate a microelectronic component 
on a printed circuit board. The top and bottom horizontal 
surfaces of the enclosure are such that their temperatures 
can be varied independently. Heat exchangers CFig. 2) used 
in conjunction with a circulating bath act as these horizon- 
tal surface. A cartridge heater housed in a stainless steel 
cover CFig. 3) acts as the 40-pin DIP model. This assembly 
is hereafter referred to as the heater. A calibrated resis- 
tor put in series with the heater and its power supply allows 
the input power to be accurately determined. Eight thermo- 
couples are placed throughout the enclosure for temperature 
measurements. Two are placed on each heat exchanger and 
four on the cartridge heater cover. 

Several criteria were established to be used as guidelines 
for design. The following is a list of these criteria and 
their implications. 



19 



HCAT EXClIANGrR 




»/i 

z 

o 



o 



%/% 

S 

U1 

z 



o 






5 h 

CD 

jQ 

B, 

03 

U 



c/^ 

CD 

H 



CD 

=5 

dO 



20 



COVER PLATE 




21 



Figure 2 . Heat Exchangers 
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Figure 3. Heater Assembly 



1 . Dimens ions 

As previously noted, the heater was to model geomet- 
rically a 40-pin DIP. Micrometer measurements of the height 
and width of a packaged 40-pin DIP were made and the rectan- 
gular heater cover is made to approximate these dimensions. 
Further, the enclosure's height was to be 6 times the verti- 
cal dimension of the heater cover and its width was to be 
between 3 and 4 times the heater cover's horizontal dimen- 
sion. This provides a good approximation to the spacing of 
a vertical printed circuit board placed in an electronic 
module . 

2 . Heater Location 

The apparatus was to be made such that the heater 
could be raised or lowered on the vertical wall and that the 
width of the enclosure could be made to vary from 1 to 2 inches. 
This would provide a wide range of enclosure sizes and heater 
locations for investigation. It was determined that the 
ability to move the vertical wall on which the heater is 
attached is better for preserving fluid-tight integrity and 
smoothness. Machine milled interconnecting parts are success- 
fully used to meet this and the variable chamber width 
requirements . 

3 . Visualization Technique 

The thymol-blue pH indicator flow visualization tech- 
nique proposed by Baker [20] was to be used. This is a non- 
intrusive technique that uses a pH indicator in distilled 
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water. When a voltage is applied across the fluid it percep- 
tibly darkens near the negative electrode due to a local 
change in pH induced by proton transfer. If there is fluid 
motion away from the cathode, as there is in this experiment, 
the darkened fluid traces the path of motion before dissipat- 
ing back to the original lighter color and neutral pH of the 
bulk fluid. In addition to its non-intrusive quality no 
density differences occur within the fluid during voltage 
application since the pH indicator remains an ion in solu- 
tion. The specifics of the fluid preparation are discussed 
in Chapter III. After much testing it was determined that 
the use of this technique requires any metal coming into 
contact with the liquid solution Ce-g*> heat exchangers, 
heater) to be relatively inactive with respect to the 
Galvanic series. This prevents corrosion as indicated by 
O’Conner [21]. To this end, a stainless steel shim strip 
running the length of the vertical wall opposite the heater 
is used as the cathode. The heater acts as the anode. 

This arrangement provides excellent visualization of the 
two dimensional flow of the darkened solution. 

4 . Other Considerations 

One eighth inch plexiglas is used as the vertical 
viewing wall to minimize optical distortion. Further, all 
vertical walls are heavily wrapped with foam rubber insula- 
tion to minimize heat loss via conduction through them 
during experimental operation. 
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B . COMPONENTS 



1 . Heater 

The heater assembly (Fig. 3) consists of an 8 inch 
long, 1/2 inch wide, 1 inch high, and 3/16 inch thick 304 
stainless steel channel that houses a narrow rectangular 
cartridge heater backed by phenolic insulation that runs 
the length of the channel. The cover has four 1/16 inch 
deep slots located on the horizontal faces 2.75 inches from 
the ends of the channel. In each of these slots a copper- 
cons tantan thermocouple is spot welded nearest the cover's 
vertical face. The thermocouple wires run out the back of 
the slots into the channel. The slots are filled with steel 
epoxy. Three evenly spaced screw holes are tapped into the 
top and bottom of the cover so that the entire assembly can 
be attached to a vertical wall in the enclosure. Prior to 
installation the cover is polished to a mirror-like finish 
using 180 to 600 grit special silicon carbide grinding paper 
and then cleaned with common rubbing alcohol. 

2 . Heat Exchangers 

The two aluminum heat exchangers measure 1.5 by 2.5 
by 11.0 inches. The top heat exchanger is seen in Fig. 2. 
They each have three rectangular channels separated by 1/8 
inch walls. The channels are machine milled into a 1.375 by 
2.5 by 9.25 inches solid piece of stock aluminum and then 
covered with a 1/8 inch thick plate secured by ten no. 4 
screws. Slots 1/8 inches deep are milled 8.0 inches apart 
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on the side opposite the cover plate to accept vertical walls. 
One slot is 3/16 inches in width and the other 3/4 inches in 
width. Both slots traverse the width of the heat exchanger. 

In addition, two 1/16 inch deep slots are milled 2.75 inches 
from the outer slots and they run to within 1 inch of the 
other side of the exchanger. These last two slots have 
copper-constantan thermocouples spot welded into them and 
are covered with steel epoxy. A 1.5 by 8.0 inch 0.001 gage 
stainless steel shim facing is epoxied in between the two 
outer slots on the side of the exchanger opposite the cover 
plate. It is positioned flush against the corner of the 
heat exchanger that is one inch from the thermocouple slots. 
This facing is the relatively inactive portion of the ex- 
changer that acts as the horizontal wall of the chamber. 

The 1.5 by 2.5 by 0.875 inch headers are also machined 
from single pieces of aluminum stock. They each have a 1/4 
inch od, 1/8 inch id exterior pipe for coupling with a cir- 
culating bath and a 2 square inch hollow interior positioned 
so as to get maximum flow in each exchanger channel. The 
headers are held to the main body of the exchangers by a 
single no. 6 screw. All seams are fitted with epoxy prior 
to assembly. 

Both heat exchangers were tested for temperature 
distribution and leakage prior to use in the experiment. 

They were coupled to a circulating bath and temperature 
measurements were taken using the two installed thermocouples 
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and an additional similar thermocouple that was randomly 
moved about the exchanger surface. The thermocouples were 
run to a digital pyrometer with an accuracy of 1 degree 
Fahrenheit. It was determined from this test that the heat 
exchangers are, upon reaching steady state, at a constant 
temperature throughout; i.e., the pyrometer reading was the 
same for each thermocouple. 

3 . Test Chamber 

The test chamber seen in Fig. 1 forms a rectangular 
enclosure measuring 1.5 inches in width, 8.0 inches in depth 
and 6.0 inches in height. The top and bottom horizontal sur 
faces of the enclosure consist of the stainless steel shim 
facing of the heat exchangers. Three of the vertical walls 
are entirely plexiglas. The wall holding the heater is 10.5 
inches wide, 24 inches high, and 11/16 inches thick. The 
heater is secured to the center of this wall with six no. 4 
screws and rubber o-rings that are counter sunk into the 
wall. In addition, a 1/2 inch diameter hole is located 1 
inch from one end of the wall to accommodate the cartridge 
heater power leads and thermocouple wires from the heater 
cover. The wall has two 1/4 inch deep slots, one 3/16 
inches in width and the other 3/4 inches in width, 8 inches 
apart and 14 inches long running parallel along its height. 
These slots receive the adjacent vertical walls. They also 
allow the wall and subsequently the heater to be raised or 
lowered while maintaining fluid tight integrity. 
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The other two plexiglas vertical walls are the ends 
that comprise the 1.5 by 6.0 inch dimensions of the enclosure. 
They are both 2 inches by 6.25 inches but one is 11/16 inches 
thick while the other is 1/8 inches thick. The thinner wall 
is used for visualization and photography and was chosen to 
be thin in order to minimize optical distortion. Both walls 
slide into the 1/4 inch slots of the two larger adjacent walls 
and into the 1/8 inch deep slots on the heat exchangers. 

The final vertical wall of the enclosure is a 10.0 
inch by 6.0 inch plexiglas frame that house a glass plate. 

The plate is installed with epoxy and will be used for future 
studies in conjunction with laser doppler velocimetry. Like 
the other large wall it has two 1/8 inch deep slots that 
receive the adjacent walls. The wall has lips that slide 
into the heat exchanger slots as well as two 1/16 inch id 
standpipes made of 1/8 inch stainless steel tubing. One 
standpipe is located near an upper corner and the other is 
diagonally opposite. They are used for filling/emptying the 
enclosure and relieving internal pressure during experimen- 
tal operation. 

C. ASSEMBLY 

Once the heater has been attached to the large wall the 
apparatus is assembled with the exception of the glass wall 
framed by plexiglas. Dow Corning 732 RTV adhesive/sealant 
is applied externally to all seams not involved with visuali- 
zation, to the power lead/thermocouple hole, and to the heads 
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of all heater attaching screws. A thin coat of Dow Corning 
3140 RTV coating and silastic is used on all internal seams 
and on all of the external seams of the smaller walls. The 
glass wall framed by plexiglas is then fitted into place 
along with a taut 1/4 inch wide strip of 0.001 gage 302 
stainless steel shim. The strip is centered internally on 
the wall and runs its height. A few inches of the strip 
traverse the seams between the wall and the heat exchangers 
(both top and bottom) to the outside of the test chamber. 

This strip acts as the cathode for the thymol-blue visuali- 
zation technique. The sealants are then applied to all 
remaining seams and the apparatus is tightened by specially 
made clamps. Fig. 4 is a photograph of the assembled apparatus. 

The vertical wall of the apparatus that holds the heater 
has to be removed, placed at the desired heater height, and 
then resealed between test runs. If removing the wall between 
test runs is not desired the apparatus can be turned upside- 
down to attain a different heater height. 

D. INSTRUMENTATION 

1 . Power to the Heater 

The heater is run in series with a precision resistor 
that was measured to have a 2.0317 ohms resistance on a 
Rosemount commutating galvonometer bridge. This arrangement 
allows a voltage measurement across the resistor to be made 
and, using its known resistance, the current through the 
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Figure 4. Photograph o£ Experimental Apparatus 
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circuit can be calculated. A measurement of the voltage 
across the heater is also taken. The product of the current 
and the voltage across the heater gives input power to the 
heater. Both voltages are measured using a Dana digital 
multimeter . 

Because of the 180 watt, 120 volt rating of the car- 
tridge heater two different power supplies are used in order 
to attain certain voltage levels. A Lambda regulated DC 
power supply is used for voltages up to 60 volts. A Variac 
voltage AC transformer is used for voltages greater than 60 
volts. A concern about using the latter was the fluctuating 
voltage provided by house current and its possible effect on 
measurements. As it happened, temperature measurements re- 
mained stable. However, voltage measurements across the 
heater fluctuated up to 3.0%. Subsequently, calculations 
of input power to the heater are less accurate with this 
mode of operation than the power calculations made when the 
Lambda regulator is used. A better method when using the 
Variac would have been to put it in series with a line volt- 
age regulator in order to eliminate voltage fluctuations 
caused by house current. 

2 . Voltage Across the Fluid 

The voltage source for the electrochemical technique 
is also a Lambda regulated DC power supply. Voltages across 
the electrodes range from 3.0 to 8.0 volts depending upon 
the contrast of the fluid required for photography. The 
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voltage is measured with the Dana digital multimeter. The 
stainless steel shim strip running the height of the verti- 
cal glass wall framed by plexiglas acts as the cathode while 
the stainless steel heater cover acts as the anode. Elec- 
trical access to the heater cover is via a wire spot welded 
to the bottom center stainless steel attaching screw that 
holds the heater assembly to the vertical wall opposite the 
cathode . 

Originally the heater was utilized as the anode and 
the top heat exchanger acted as the cathode. This arrange- 
ment proved unacceptable because flow paths could not easily 
be discerned from the viewing perspective. The inclusion of 
the shim strip as the cathode solved this problem by limit- 
ing the darkened flow to a single path located in the center 
of the enclosure. 

3 . Temperature Measurement 

Eight, 30 gage copper-constantan (teflon/teflon in- 
sulation) thermocouples dispersed throughout the enclosure 
are used for temperature measurements. Two thermocouples 
are located on each horizontal face of the heater cover and 
two on each heat exchanger. All thermocouples were fabricated 
using a Dynatech welder and calibrated on a Rosemount commu- 
tating bridge in conjunction with a platinum resistance tem- 
perature standard and calibration bath. 

The thermocouples run to a Newport digital pyrometer 
model 267 where temperature readings are made. The pyrometer 
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gives readings to the nearest degree Fahrenheit. If, during 
the measuring process, a reading fluctuates between two 
temperature extremes it is assumed to be the average of 
the two. Finally, the ambient temperature is measured with- 
in two feet of the apparatus with another thermocouple run 
to the pyrometer. 
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III. EXPERIMENTAL PROCEDURE 



A. ELECTROCHEMICAL FLUID PREPARATION 

As noted earlier, the thymol-blue pH indicator flow vis- 
ualization technique as described by Baker [20] is used. 
Further, O'Conner suggests using thymol-blue sodium salt 
vice pure thymol-blue Cthymolsulphonephthalein] crystals in 
the solution to avoid filtering of undissolved crystals [21] . 
This advice is followed in the preparation of the fluid. 
First, 0.2 grams of thymol-blue sodium salt are dissolved 
in 1.0 liter of distilled water. This gives a 0 . 02 % by 
weight light colored solution with no undissolved solids. 
Next, 0.1 normal sodium hydroxide is added drop by drop to 
the solution until it titrates to a dark blue color. The 
dark color is an indication of a slightly basic pH. Finally, 
the solution is titrated back to its original lighter color 
by adding drop by drop 0.1 normal hydrochloric acid. The 
fluid is now slightly acidic and very close to its neutral 
pH value of 7. After preparation, a small amount of fluid 
is poured into a shallow glass container. Wires attached to 
a 1.5 to 7.5 volt dry cell battery are used as the electrodes 
and the fluid is tested for local color change. All tests 
indicate adequate color change. 
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B. APPARATUS PREPARATION 



After the sealants are cured, 1/4 inch tygon tubing is 
attached to the heat exchangers using stainless steel hose 
clamps. The exchangers are then interconnected in series to 
a refrigerated circulating bath placed six feet from the 
apparatus. The top heat exchanger is placed upstream of the 
bottom exchanger. Common tap water is used as the refriger- 
ated circulating coolant. The apparatus is then laid hori- 
zontally on its largest wall. 1.5 feet of 1/8 inch tygon 
tubing is attached to each standpipe and the enclosure is 
filled with the thymol-blue solution via a simple siphon 
system. Manual manipulation of the apparatus and standpipe 
tubing allows most bubbles to be driven out of the enclosure. 

A specially made steel leveling stand with four height 
adjustable legs and a slot for the largest vertical wall had 
been made to hold the apparatus. The stand is placed on a 
piece of plywood under which lies two inches of foam rubber 
insulation. This system serves to damp out any vibrations 
occurring in the laboratory. The apparatus is placed verti- 
cally upright onto the leveling stand and secured in place. 

The two larger vertical walls are then layered with 3/4 inch 
thick insulation. After insuring that the view through the 
test chamber to a white background is unobstructed, 3/4 inch 
thick insulation is placed on the remaining vertical walls. 

The apparatus is then leveled both horizontally and vertically 
by means of a bubble level placed on the top heat exchanger 
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in two orthogonal directions and by concurrent adjustment of 
the stand legs. This leveling process is crucial for ensur- 
ing horizontal and vertical alignment of the chamber. Fig. 5 
is a photograph of the assembled apparatus and instrumentation. 

C. TEST PROCEDURE 

1 . Heater Location 

Test runs were made for three heater locations: (1) 

the top horizontal face of the heater located 1.0 inch from 
the top of the enclosure, (2) the heater centered on the 
vertical wall of the enclosure, and (3) the bottom horizontal 
face of the heater located 1.0 inch from the bottom of the 
enclosure . 

2 . Initial Instrument Settings 

Temperature measurements of all thermocouples are 
made to insure that the apparatus and its liquid contents 
are at a single temperature. If this criteria is met the 
circulating bath to the heat exchangers is turned on to a 
predetermined temperature. Shortly thereafter the power 
supply to the heater is turned on and the voltage slowly 
increased to another predetermined value that results in 
the approximate desired temperature difference between the 
exchangers and the heater. The values of the settings for 
the circulating bath temperature and the voltage to the 
heater are determined soley from experience garnered from 
a multitude of test runs . 
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Figure 5. Photograph of Experimental Apparatus and Instrumentation 





3 . Instrument Readings 

Thermocouple temperature measurements and voltage 
measurements across the heater and resistor are taken every 
twenty minutes until successive temperature measurements 
remain unchanged to within the accuracy o£ the pyrometer. 
Measurements are then taken every ten minutes. It is assumed 
that steady state is attained within the enclosure i£ no 
thermocouple measurement changes by more than one degree 
during any ten minute period. This process takes between 
two to £our hours depending on the initial temperature di£- 
£erence between the heater and the exchangers. The smaller 
temperature di££erence takes the least amount o£ time. 

Once steady state is attained the power supply that 
creates the voltage across the electrochemical £luid is 
turned on. Using a digital multimeter as an indicator, the 
voltage is increased £rom 0.0 to 3.0 volts. The insulation 
is then removed £rom the two smaller vertical walls and the 
£luid voltage is varied until a visually suitable path o£ 
darkened £luid ensues. 

4 . Photographic Techniques 

A Nikon Ftn with a 50 mm £2 lens is used £or all black 
and white photography. A Canon Ftb with a 50 mm £1.4 lens 
is used £or all color photography. Black and white photo- 
graphs are taken with Kodak Tri-x pan asa 400 £ilm; color 
photographs are taken with Kodak Ektachrome asa 64 daylight 
£ilm. Black and white photographs are taken £irst. Both 
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cameras are secured to a tripod. The Nikon is placed between 
14 and 15 inches from the 1/8 inch thick plexiglas vertical 
wall. A combination of f-stop settings of 8 and 11 coupled 
with shutter speeds of 1/60 and 1/125 of a second are used. 

The Canon is placed at a distance of between 12 and 13 inches 
from the viewing wall. Combinations of f-stop 2 and 4 with 
shutter speeds of 1/4 and 1/15 of a second are used with 
the Canon. 

Two L-shaped pieces of cardboard are attached together 
so as to create a rectangular opening approximately the size 
of the small vertical wall. This is placed between the camera 
and the enclosure. A cloth is hung from this cardboard aper- 
ture and draped over the camera during photography to elimi- 
nate camera reflection on the plexiglas. Another larger 
piece of cardboard covered with white paper is placed 1 foot 
from the rear small vertical wall to act as a backdrop. 

Finally a 650 watt colortran multi-6 set lamp is positioned 
2 feet from the white cardboard backdrop and 1.5 feet from 
the apparatus. The lamp is directed towards the white back- 
drop and turned on only when photographing the flow. 

A major concern was that the heat generated by the 
lamp would disturb the flow. However, due to the insulation 
remaining between the lamp and the enclosure and the short 
interval of time the lamp is on the viewed flow remains two 
dimensional. This scheme provides an excellent view down the 
length of the enclosure to a white background that adequately 
contrasts the darkened fluid. 
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D. DATA ANALYSIS 



There are six length dimensions within the enclosure, any 
one of which could be chosen as the characteristic length, 

L, for determining the Nusselt number and Rayleigh numbers. 

The 1.0 inch vertical face of the heater is chosen since it 
is the surface that not only adds the majority of energy to 
the system but it also is oriented in the direction of the 
gravitational driving force of free convection. The proper- 
ties of the fluid are evaluated at the average of the hot 
surface (heater) and cold surface (heat exchangers) temperatures. 

1 . Determination of Nusselt Number 

The Nusselt number is defined as follows; 

Nu = (1) 

Kf 



To evaluate the heat transfer coefficient an estimate has to 
be made of the heat lost through the vertical walls of the 
enclosure. A parallel thermal resistance circuit for con- 
duction through the three vertical walls in contact with the 
heater is employed for this purpose. Heat loss through the 
vertical wall opposite the heater, all contact resistances, 
and conduction via the heater assembly attaching screws are 
assumed negligible. Knowing the total thermal resistance 
for conduction the heat loss is calculated using a variation 
of Fourier's Law of Conduction: 



Q 



cond 




^cond 



(2a) 
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where 



R 



^ L 
cond ^ 



C2b) 



The length of the conduction path, L, for all materials is 
known. The thermal conductivity, k, for each material is 
determined using a table of properties [22], The area. A, 
is assumed to be that of the cartridge heater perpendicular 
to the conduction path. Finally, the temperature difference 
is assumed to be the difference between the ambient, T„ , 
and the heater cover temperature, T^ ; this further neglects 
the convective thermal resistance on the outside of the in- 
sulation. In the event T is less than T , heat loss is 
assumed to be zero. 

After calculating the energy lost via conduction and 
having already determined the energy into the system 
from the voltage measurements across the heater and the pre- 
cision resistor a simple energy balance is used to determine 

the energy convected into the fluid CQ ) : 

=>•' ^conv 



Q 



conv 




C3a) 



where 




(3b) 



and 



= voltage across heater 
V2 = voltage across resistor 



R = 2.0317 Ohms 
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Then, using Newton's Law of Cooling, 



Q 



conv 



h A (T^ - 



C4) 



h is calculated. ?£ is chosen to be the temperature of the 
heat exchangers which are the ultimate heat sink for the 
enclosure. The exterior surface area of the heater cover 
is A. With all variables in equation (1) now known the 
Nusselt number is calculated. 

To ensure the accuracy of this method of analysis a single 
iteration is performed on the Nusselt number by use of another 
thermal resistance for conduction circuit. This time a single 
circuit with a lone resistance is run through the steel heater 
cover to the cartridge heater with the intent of determining 
the cartridge heater temperature. Using the heater cover 
temperature (T^) , length of the conduction path through the 
heater cover CL) , the thermal conductivity for 304 stainless 
steel Ckgg) , the exterior area of the heater cover CA) , and 
the previously calculated convection energy into the enclo- 
sure (Q ) , it is a simple matter to determine the cartridge 
^conv ^ ° 

heater temperature CTj^) from a rearrangement of Fourier's Law 
of Conduction: 



L Q 




( 5 ) 
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Tg in equation (2a) is then set equal to Tj^ and a single 
iteration is performed to re-evaluate the Nusselt number. 
Iterative analysis of Nusselt numbers for all test runs 
prove to fall within 3.0% of the initial calculation thus 
eliminating the necessity for further iterations. A sample 
calculation is contained in Appendix A, 

2 . Determination of Rayleigh Number 

One definition of the Rayleigh number is as follows: 



where 



Ra 



Gr 



(Gr) (Pr) 

g B (Tg - T^) 

2 

V 



Pr = 



V 

a 



(6a) 

(6b) 

(6c) 



In the calculation of the Grashof number (Gr) , the same 
characteristic length (L) is used as that used to determine 
the Nusselt number; that is, the 1.0 inch vertical face of 
the heater cover. In addition, the same temperature differ- 
ence is used as that in equation (4) (i.e., the average 
heater cover temperature minus the average heat exchanger 
temperature) . The remaining terms in equation (6b) are eval 
uated at the average temperature of all thermocouple measure 
ments by a cubic polynomial that had been fit to data from 
a table of properties [23] for temperatures of water ranging 
from 50 to 100 degrees Fahrenheit. 
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The Prandtl number (Pr) is also evaluated by a cubic 
polynomial that had been fit to data from the same table 
of properties, at the same temperature, and for the same 
working fluid. The Rayleigh number is then calculated. 

A sample calculation is contained in Appendix A. A sample 
uncertainty calculation is contained in Appendix B. 
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IV. RESULTS 



A. VISUAL 

Photographs of the natural convection flow within the 
enclosure for a representative number of heater locations 
and Rayleigh numbers are presented in Figs. 6, 7, and 8. 

Figure 6 depicts flows with the top horizontal face of the 
heater located 1.0 inches from the top of the enclosure. 

Figure 7 depicts flows with the heater centered in the en- 
closure. Finally, Figure 8 depicts flows with the bottom 
horizontal face of the heater located 1.0 inches from the 
bottom of the enclosure. 

Some general observations can be made about the natural 
convection flow within the enclosure. First, the flow is 
dual-celled. There is a buoyancy driven upper cell and a 
shear driven lower cell. The location of the cell separa- 
tion layer is dependent upon the heater location within the 
enclosure; as the heater is raised so, too, is the layer. 

Secondly, the flow of the buoyant driven upper cell fol- 
lows the geometry of the enclosure. The fluid comes off 
the vertical face of the heater, up the vertical wall to 
which the heater is attached, across the cold heat exchangers, 
then down the vertical wall opposite the heater. Finally, 
the fluid moves along the cell separation layer where it 
interacts viscously with the lower cell and back towards the 
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Figure 6* Photographs of Flow with Heater Near Top of Enclosure 
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Figure 7. Photographs of Flow with Heater Centered in Enclosure 









t 
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Figure 8. Photographs of Flow with Heater Near Bottom of Enclosure 





heater. No definitive comments can be made about the total 
motion of the fluid in the lower cell; only a very slow moving 
flow was observed slightly below the separation layer. This 
flow was crossing the width of the enclosure towards the ver- 
tical wall holding the heater. 

The photographs also show flow being generated into each 
cell. This is the result of a recurring problem encountered 
using the Baker flow visualization technique. Baker suggests 
applying a maximum of 5.0 volts across the fluid when elec- 
trode spacing is 0.5 cm. These limitations preclude genera- 
tion of bubbles near the cathode. Utilizing the limits imposed 
and a network of cathodes the ensuing flow can be seen as it 
moves away from the area of interest. In this experiment a 
single cathode is used. The area of interest is the entire 
enclosure. Hence, it is necessary to apply up to 8.0 volts 
in order to observe the flow throughout the entire enclosure. 
Bubbles are therefore generated that cause convective mass 
transfer up the stainless steel shim strip. The buoyant flow 
circulation of the upper cell collides with the mass transfer 
associated with the convected bubbles resulting in the intru- 
sion of flow into the cells. Efforts were made to minimize 
this effect by careful adjustment of the voltage during test 
runs . 

The distinct difference between the photographs of Fig. 7 
clearly illustrates the advantage of utilizing a stainless 
steel shim strip as the cathode for the visualization technique 
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As previously noted, the heater was originally utilized as the 
anode while a heat exchanger acted as the cathode. This pro- 
vided an excellent darkening of the fluid but also proved in- 
adequate in terms of flow visualization (tight photograph of 
Fig. 7). The inclusion of the shim strip as the cathode allows 
visualization of the two dimensional flow that occurs in the 
middle of the enclosure (left photograph of Fig. 7). 

B. QUANTITATIVE 

Table I summarizes the results of calculations using data 
obtained for each test run. Based on these calculations the 
Nusselt and Rayleigh numbers were plotted on a log- log graph 
(Fig. 9). Since the temperature measurement uncertainty in- 
herent in the experiment produces Nusselt and Rayleigh number 
uncertainties ranging from five to twenty percent (see Appendix 
II) it is difficult to draw definitive conclusions from the 
data. However, there does appear to be a trend suggesting 
that as the heater is raised within the enclosure the Nusselt 
number decreases. 
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TABLE I 



Summary of Results 



Tg (°F) 




Tf (°F) 


Qin (Btu/hr) 


Nu 


Ra (x 


Heater 


Near 


Top 








65.75 




54.00 


26.41 


4.95 


0.99 


70.25 




50.00 


55.75 


6.12 


1.73 


72.25 




38.88 


102.68 


6.84 


2.17 


81.25 




38.88 


146.85 


7 . 65 


3.61 


82.50 




39.25 


145.52 


7 . 38 


3.85 


89.75 




36.75 


197.58* 


8.19 


5.33 


Heater 


Centered 








66.75 




57.00 


18.94 


4.23 


0.92 


69.25 




59.00 


25.99 


5.58 


1.07 


68.50 




57.50 


28.79 


5 .76 


1.09 


69.00 




50.50 


51.55 


6.12 


1.55 


66.75 




42.00 


68.15 


6.12 


1.49 


70.50 




40.63 


99.01 


7.38 


1.94 


70.50 




39.25 


103. 80 


7.38 


1.95 


77.75 




39.75 


144.89 


8.46 


3.00 


75.50 




36.75 


147.37 


8.46 


2 . 61 


79.00 




39.50 


147.33 


8.28 


3.21 


80.00 




39.25 


148.06 


8.01 


3. 38 


88.00 




40.00 


185.42* 


8.46 


5.01 


99.50 




38.00 


248.91* 


8.73 


7.95 


Heater 


Near 


Bottom 








68.25 




56.50 


28.78 


5.40 


1.13 


68.00 




51.00 


51.24 


6.66 


1.40 


68.50 




40.00 


98. 12 


7.74 


1.70 


75.75 




40.00 


139.54 


8.64 


2.68 


Note : 


(1) 


* indicates 


Variac used as 


power supply 
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Figure 9. Data Plot of Nusselt versus Rayleigh Number 



V. CONCLUSIONS 



The author has found no previous natural convection study 
of a liquid filled rectangular enclosure with a heated pro- 
trusion attached to a vertical wall, temperature variant 
horizontal walls, and the remaining walls adiabatic. However, 
some investigators using similar configurations have reached 
the same conclusions about to be presented. 

A. VISUAL 

Two cells occur within the enclosure during the natural 
convection flow. The upper cell is driven by buoyant forces 
and essentially follows the geometry of the enclosure. The 
lower cell is shear driven. The location of the cell separa- 
tion layer is dependent upon the height of the heater in the 
enclosure . 

B. QUANTITATIVE 

Calculations indicate a trend for the Nusselt number to 
decrease as the height of the heater within the enclosure is 
increased. It is interesting to note that for the Rayleigh 
numbers calculated in this experiment (on the order of 10^ 
to 10^) Chu, et al [11] plot this same tendency for an air 
filled square enclosure that has a heater imbedded in a ver- 
tical wall and cooled horizontal walls. It is Chu's inves- 
tigation that most closely resembles this experiment. 



53 



VI. RECOMMENDATIONS 



A. IMPROVEMENT TO EXPERIMENT 

1 . Apparatus 

Fluid leakage into the heater assembly occurred during 
the course of the experiment. This caused the cartridge heater 
to fail. Replacement of the cartridge heater necessitated dis- 
mantling the apparatus. To avoid the problem of cartridge 
heater malfunction, either a liquid immersible cartridge heater 
should be used or else an improved design of the heater cover 
should be sought. 

Though the photographs proved satisfactory it is fur- 
ther suggested that experimentation with various camera filters 
be done. The contrast of the observed flow may improve to the 
point where excessive voltage across the fluid may not be re- 
quired, thereby, eliminating the bubble generation effect. 

2 . Instrumentation 

In order to lower the experimental uncertainty it is 
suggested that a more accurate temperature measuring instru- 
ment be used for thermocouple readings. Further, to preclude 
the use of the Variac transformer a cartridge heater with the 
same power (180 watts) but rated at 60 instead of 120 volts 
should be used. This will allow all test runs to be made 
with the Lambda power supply as the power supply to the 
heater . 
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B. ADDITIONAL EXPERIMENTAL WORK 



It is suggested that the following areas of study be 
experimentally explored: 

1. Decrease the enclosure width to see if any additional 
cells develop within the enclosure. 

2. Observe the effects of independently varying the heat 
exchanger temperatures. In conjunction with this, do not 
allow the coolant temperature to go below fifty degrees Fah- 
renheit to avoid the density inversion that occurs in water 
at forty degrees. As can be seen in Table I, the heat ex- 
changer temperatures for some test runs fell below this value. 

3. Continue to use the stainless steel shim strip as the 
cathode but seek methods of decreasing the effect of the 
bubbles generated when voltage is applied across the fluid. 
Alternatively, design a network of cathodes (as suggested by 
Baker) that encircles the enclosure. 
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APPENDIX A: SAMPLE CALCULATIONS 



A. DETERMINATION OF INPUT POWER 

The dimensions of the cartridge heater are 8.0 by 0.25 
by 0.725 inches. Table II is the data that was recorded 
for a representative run on 3 November 1983 with the heater 
located near the bottom of the enclosure. Refer to Fig. 1 
for thermocouple locations. 

TABLE II 

Test Run Data for 3 November 1983 

40 40 40 40 73 77 77 76 71 58.51 1.42 

Notes: (1) all temperatures in degrees Fahrenheit 

C2) all voltages in volts 

Using equation 3b the input power is calculated to be: 

Q. = (3.4123 Btu/hr-Watt) = 139.54 Btu/hr 

(2.0317 ohms) 

B. NUSSELT NUMBER 

The heat loss via conduction through the vertical walls 
is determined by constructing a thermal resistance network 
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(Fig. 10). The thermal conductivity of each material 
(Table III) is determined from a table of properties [22], 



RUBBER 

INSULATION PLEXIGLAS 




Figure 10. Thermal Resistance Network for Calculation of Losses 
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TABLE III 



Thermal Conductivity of Materials 



Material 

rubber insulation 
plexiglas 
phenolic 
304 steel 



K (Btu/hr-ft-°F) 
0.0225 
0.0821 
0.0175 
8.8693 



Table IV is a summary of the individual resistances for Fig. 
10 that were determined using equation 2b. 



TABLE IV 



Resistance 




Thermal 

Path(l) 


Resistances to Conduction 
Material Value 


(°F-hr/Btu) 






a 


plexiglas 


116.9 


^2 




a 


rubber insulation 


2562.7 


R, 




b 


plexiglas 


643.2 


R4 




b 


rubber insulation 


2562.7 


‘'s 




c 


phenolic 


8.7 






c 


plexiglas 


20.1 


R? 




c 


rubber insulation 


80.1 


Note: (1) 


conduction path through 

a. viewing wall 

b. wall opposite viewing 

c. wall to which heater is attached 
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The equivalent 



resistance for 

R = 2679.6 
a 

= 3205.9 

R = 108.8 

c 



each path is: 
°F-hr/Btu 
°F-hr/Btu 
°F-hr/Btu 



Therefore , 



R 



cond 



(R^) (Rt,) (Rc> 

(R^)CRi,)*(R3)(R^)*CRb)(R,) 



101.28 °F-hr/Btu 



Employing equation 2a, 



Q 



cond 



(75.75 - 40.0) 
(101.28) 



0.047 Btu/hr 



The Nusselt number is determined by first using equation 
3a to calculate the energy convected into the fluid: 

Qconv " (139.54 - 0.047) = 139.495 Btu/hr 

Then, the heat transfer coefficient is calculated as (equa- 
tion 4) 

^ , (139.495) (144) , 35 

(75.75-40.0) ( 16 ) 

The thermal conductivity of water is determined from a table 
of properties [23] , 
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0.3383 



k 



£ 



Btu/hr-ft-°F 



Finally, the Nusselt number can be calculated (equation 1 ) , 



Nu = = 8,64 

(0.3383) 



An iterative check is performed by determining the car- 
tridge heater temperature (equation 5) as 

T ^ (3/192) (139.495) ^ 75 75 = 77.96 °F 
^ (8.659) (1.33) 



Reapplication of equations 2a, 3a, 4, and 1 respectively with 
set equal to Tj^ yields 



Q 



cond 



(77.96 - 71.0) 
(101.285) 



0,0687 Btu/hr 



Qconv " (139.54 - 0.0687) = 139.474 Btu/hr 

h = (139.474) (144) — ^ 35,10 Btu/hr-ft^-°F 

(77.96 - 71.0) (16) 



From the table of properties [23] the fluid conductivity 
remains 



k^ = 0.3383 Btu/hr/£t/°F 



60 



Therefore , 

Nu - g /m . 8.64 

(0.3383) 

The percent difference in Nusselt numbers after one itera- 
tion is zero percent. 

C. RAYLEIGH NUMBER 

Using the property table [23] at the average temperature 
of the heater and heat exchangers the Prandtl number and 
Grashof number can be calculated as ; 

Pr = 8.322 

^ = 15.58 X 10^ l/°F/cubic ft 
Gr = (15.58) (77.96 - 40.0) (1/1728) = 322.4 x 10^ 

Therefore, the Rayleigh number can be determined using 
equation 6a: 

Ra = (322.4 x 10^) (8.322) = 2.683 x 10^ 
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APPENDIX B: UNCERTAINTY ANALYSIS 



Uncertainties in the experiment are evaluated by employing 
the root mean square method of Kline and iMcClintock [24]. 
Table V is a summary of measurement and property values and 
their associated uncertainties: 



TABLE V 

Uncertainty in Variables 



Variable 


Uncertainty 


Basis for Uncertainty 


g6/v^ 


1/250 


assume last digit 


kf 


1/300 


assume last digit 


k 

g 


5.0 1 


[22] 


k. 

1 


7.0 1 


[22] 


k 

P 


7.0 % 


[22] 


k 

ss 


5.0 1 


[22] 


L (1) 


0.001 in 


resolution of measurement device 


Pr 


1/600 


assume last digit 


R 


0.0001 ohms 


calibration on galvonometer 


T (2) 


1.0 F 


resolution of measurement device 


V^(3) 


0.01 volts 


resolution of measurement device 


V, (4) 


2.0 volts 


resolution of measurement device 




0.01 volts 


resolution of measurement device 



Notes: (1) all dimension measurements 

(2) all thermocouple temperature measurements 

(3) Lambda regulator as power supply 

(4) Variac as power supply 
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The following variable uncertainties are determined using 
the aforementioned RMS equation on the data presented in 
Table II of Appendix A. 



WQin = [ ( 



(WV,)(V^) (WV2)(V, ) (WR)(V,HV2) " 

—) + ( ^ —) + C —) ] 



1/2 



R 



R 



R 



= ^ ^ (O.OlHl.42) ^^, ^- (0.01) (58.51) ^^, ^ (0.0001) (58. 51) (1.42)^ 
(2.0317) (2.0317) (2.0317)^ 



V2 



= 0.983 Btu/hr 



WQ 



cond 



WT 2 
[(— ) + ( 



R 



a ^ ^ (. cond _a s_^ ^ 



1/2 



cond 



R 



cond 



^^cond^‘ 



= ^ (2) (l/2) j^ ^ ( (2) (1/2) ^^ ^ ^ (0.00005) (2) ^^^ 
(101.28) (101.28) (101.28)^ 



1/2 



= 6.017 Btu/hr 



WQ = [(WQ j) + (WQ. ) ] 
^conv ^ ^cond ^ ^in"^ ^ 



1/2 



= [(0.017)^ + (0.983)^] 



1/2 



= 0.983 Btu/hr 
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1/2 



WNu = + (( WL)Ch) ^^ ^ ,( Wkf)(h)(L) ,^j 

Ck.) 



Ck^) 



Ck^)2 



^^ (2.95) (0.085) ^“^ ^ CO. 001) (35 . 10) ^ ^ , ^- C0.0QlK55.10) (0.085) ^^ j 
(0.338) (0.338) (0.338)^ 



V2 



= 0.727 



WRa = [((WGr) (Pr))^ + ( (WPr) (Gr ) ) ^ ] 



1/2 



= [((21.0 X 10-^) (8.322))^ + ((0.0)(322.0 x lO^^)) ] 



1/2 



= 178.85 X 10' 



The above uncertainty calculations are intended to be 
representative of the method used; calculations performed on 
every variable are not included. A summary of quantitative 
results for Nusselt and Rayleigh numbers and their respective 
uncertainties in relation to heater temperature minus ex- 
changer temperature for all test runs is given in Table VI. 
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TABLE VI 



Uncertainty in Experiment 



T -T. 
s f 


Nu 


Uncertainty 


Ra (x 10^) 


% Uncertainty 


Heater 


Near Top 








11.75 


4.95 


18.20 


0.99 


17.40 


20.25 


6.12 


11 . 73 


1.78 


10.52 


33.38 


6.84 


8.69 


2.17 


7 . 00 


42.38 


7.65 


7 .85 


3.61 


5.95 


43.25 


7.38 


7.79 


3.85 


5.88 


53.00* 


8.19 


7 . 89 


5.33 


5.23 


Heater 


Centered 








9.75 


4.23 


21.57 


0.92 


20.83 


10.25 


5.58 


20.55 


1.07 


19.85 


11.00 


5.76 


19.28 


1 . 09 


18 . 53 


18 . 50 


6.12 


12.53 


1 . 55 


11.40 


24.75 


6.12 


10.26 


1.49 


8.86 


29.88 


7.38 


9.19 


1.94 


7.61 


31.25 


7.38 


8.97 


1.95 


7.35 


38.00 


8.46 


8.19 


3.00 


6.39 


38.75 


8.46 


8.13 


2.61 


6 .31 


39.50 


8.28 


8 . 07 


3.21 


6.23 


40.75 


8.01 


7.97 


3.38 


6.10 


48.00* 


8.46 


8.07 


5.01 


5.52 


61.50* 


8.73 


7 . 53 


7.95 


4.87 


Heater 


Near Bottom 








11.75 


5 . 40 


18.19 


1.13 


17.40 


17.00 


6.66 


13.37 


1.40 


12.31 


28 . 50 


7 . 74 


9.43 


1.70 


7 . 90 


35 .75 


8.64 


8.41 


2.68 


6.67 



Note: (1) * Indicates Variac used as power supply. 
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